Resolving distribution of relaxation times in Poly (propylene glycol) on the crossover 

region 
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In this paper, a recently developed numerical technique [Tuncer E and Gubanski S M, IEEE 
Trans Did El Insul 8(3)(2001) 310-320] is applied to poly (propylene glycol) complex dielectric 
data to extract more information about the molecular relaxation processes. The method is based 
on a constrained-least-squares (c-lsq) data fitting procedure together with the Monte Carlo (mc) 
method. We preselect the number of relaxation times with no a-priori physical assumption, and use 
the Debye single relaxation as "kernel", then the obtained weighting factors at each MC step from 
the C-LSQ method builds up a relaxation time spectrum. When the analysis is repeated for data 
at different temperatures a relaxation-image is created. The obtained relaxation are analyzed using 
the Lorentz (Cauchy) distribution, which is a special form of the Levy statistics. In the present 
report the (3 and a relaxations are resolved for the PPG. A comparison of the relaxations to those 
earlier reported in the literature indicate that the presented method provides additional information 
compared to methods based on empirical formulas. The distribution of relaxation times analysis is 
especially useful to probe the crossover region where the a- and /3- relaxations merge and the results 
show that the relaxation after the crossover region at higher temperatures is Arrhenius-type as the 
/3-relaxation. Moreover, this relaxation is more likely to be the continuation of the /3-relaxation, but 
with a different activation energy. 

PACS numbers: 02.30.Zz, 02.70.Wz, 64.70.Pf, 77.22.-d, 77.22. Gm, 78. 55. Or 
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I. INTRODUCTION 

Traditional methods for analyzing dielectric data as 
well as impedance data consist of applying curve-fitting- 
algorithms using empirical formulas [see for example 
Macdonald^]. By such analysis information on conduc- 
tivity, molecular relaxation, liquid-glass transition, etc., 
of materials are obtained which are challenging subjects 
in condensed matter physics. A possible new overture 
that has previously been avoided by researchers and in- 
strument manufacturers has been the distribution of re- 
laxation times (drt) approach. Recently a computation 
method has been developed by Tuncer and Gubanski^ 
that can resolve unique drt for dielectric data, see for ex- 
ample Tuncer and Gubanski^, Tuncer et al.'^ and Tuncer 
et al.^. The method is based on the Monte Carlo (mc) 
technique^ and the constrained-linear-least-squares (c- 
lsq) algorithm^. The procedure obeys the Kramer- 
Kronig relations'". It is, therefore, distinct to acquire 
ohmic-losses even if their contributions are not visible 
to the naked eye. Finally, once the DRT is known it is 
trivial to obtain the time-domain relaxation function-i/ie 
responce function. 

None of the conventional techniques can resolve dipo- 
lar dielectric relaxations from measured data that con- 
tain contributions of several processes and because of 
measurement limitations. So in this paper, the devel- 
oped numerical method^ is applied to analyze dielectric 
data of atactic poly(propylene glycol) (ppg) (see Fig. 1). 
This polymer has a branched structure which is differ- 
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FIG. 1: Chemical structure of poly(propylene glycol) with the 
average molecular weight of 4000, the number of structural 
units are approx. 69. The asterisk (*) indicates the chiral 
carbon. In atactic PPG the configuration (R or S) of the chiral 
carbon is randomly distributed along the chain. 



ent from the more commonly used poly(ethylene gly- 
col). A random distribution of the absolute configura- 
tion of the tertiary carbon leads to a large number of 
different sequences, and therefore, to a complete lack of 
crystallinity also at the lowest temperatures. Anyhow, 
the orientation of the structural unit in the chain, quite 
regularly head to tail, leads to a cumulative dipole mo- 
ment along the chain contour^ (Stockmayer type-A poly- 
mer). In PPG, several relaxation phenomena are observ- 
able in the dielectric spectroscopy^"^"*. The relaxations 
are labeled as a', a and /3 in dielectric data. The a'- 
relaxation, also known as the normal-mode relaxation, is 
strongle dependent on molecular weight, and its ampli- 
tude is lower than the other ones. The a-relaxation can 
be related to the cooperative movements of the segments 
of the macromolecular chains and its temperature depen- 
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FIG. 2: Crossover region in the Arrhenius diagram^'^. a is 

the high temperature process, (5 the local Arrhenius process 
(Johari-Goldstein process), and a the cooperative process. 



dence. The /3-relaxation is suggested to be associated 
with locahzed motions of the structural units and it has 
a weaker Arrhenius-type temperature dependence^^'^^. 
Recent studies of dielectric relaxation in glass-forming 
liquids have shown a complex phenomenon, which is 
known as the crossover region^^'^^'^^"^^. In this region, 
a- and /3-relaxations merge in the Arrhenius plot (see 
Fig. 2), after the merging point, one relaxation survives 
whose origin is sought. The crossover region in glass 
forming liquids can in general be described as a transi- 
tion window from a molecular liquid at high t(UTipcratures 
(o process) to a cold liquid where cooperative process {a 
process) is present^^ as presented in Figure 2 where Tc 
and logWc are temperature and angular frequency at the 
crossover. In other words at high temperature, above the 
crossover region, one can observe the molecular charac- 
teristic relaxations, i.e. the local rearrangements, while 
below that region one can see two different processes; 
associated with local and collective processes. We there- 
fore focus on the crossover region in PPG, and probe this 
region with the DRT analysis in this study. 



II. BACKGROUND ON DIELECTRIC DATA 
ANALYSIS 

When dielectric data are obtained over a wide fre- 
quency range, not only contributions of relaxation po- 
larization, direct-current (dc) conduction and low fre- 
quency dispersion (lfd) but also contributions of space 
charges, electrode effects etc., may be present. If only 
the contributions of relaxation polarization, DC conduc- 
tion and LFD are taken into account, it is possible to dis- 
tinguish the DC conduction among the other processes. 
Kramers-Kronig relations'^ are often used for this inten- 
tion. However this requires that, (i) the dielectric system 
gives a linear response and {ii) the frequency range should 
be wide enough for the integration. The second condition 
is not usually possible for most experiments and equip- 
ment. There are several other ways to separate the DC 



conductivity contribution such as using pairs of Debye 
functions or Curie von Schweidler power-law functions- 
which are in principle the same as Kramer-Kronig trans- 
formations and finally performing measurements in the 
time domain^'''^^ . 

Besides, data analyses can be carried out in the 

frequency domain as well as in the time domain 
using, for example, complex non-linear-least-squares 
(cnlsq)^'^^"^^ or DRT methods^'^'^"^'^. In the former tech- 
nique, an analytical expression of the model is needed. 
The expression can be in the form of an equivalent 
circuit^ or of known relaxation functions^^''^^ or even as 
an addition of several functions'^^"^'*. However, Factor 
et al. has stated that for a narrow range of data cnlsq 
is not reliable. 

When there is no analytical model available for the 
relaxation behavior or single relaxation processes are 
of importance DRT can be preferred. Approaches with 
known-distribution-functions for DRT have been proposed 
in the past^^"^^. The most general ones are either 
Havrialiak-Negami^'' relaxation polarization or its deriva- 
tives^^; Cole-Cole^" and Cole-Davidson^^. There have 
been several publications on the physical bases of these 
analytical functions'^''*^"^'* . There have also been a cou- 
ple of critical comments in the literature regarding the 
DRT and its meaning, i.e., "...We have, however, no 
actual proof of necessity for the existence of a distribu- 
tion of relaxation times for a given system, nor have we 
any right to assume that a possible distribution should 
be Gaussian'"^^ and "...such theories as distribution of 
relaxation times which are purely phenomenological and 
are not based on any rigorous physical arguments. There 
are also empirical mathematical functions which repre- 
sent the experimentally observed results to a certain ex- 
tent, but which are not in any way related to the physical 
nature of the phenomena in question"'^^ . When more re- 
laxations coincide in energy or they overlap in part, re- 
searchers try to choose appropriate compositions, tem- 
perature windows and molecular weights to avoid the 
coincidence, but this is not always possible. The ana- 
lytical functions^°'''^'^^, on the other hand, are not often 
capable of resolving overlapping relaxations^^. The men- 
tioned functions are extensively used to obtain average 
or most probable relaxation times. As a result, the DRT 
have been more or less avoided as an established dielec- 
tric data analysis tool, though there have been couple of 
numerical attempts^®"^'''^®. A mathematical method to 
resolve the drt has also been proposed by Schiifcr et al. 
and Winterhalter et al. who used Tikholov regulariza- 
tion algorithms. CNLSQ has also been applied to resolve 
the DRT^^~^^. In addition arbitrary analytical distribu- 
tions were also proposed and used in the cnlsq algo- 
rithms^*^'^^'^^. As they mention, this ill-posed problem 
can not be solved using conventional LSQ and linear pro- 
gramming but in this present paper we show the useful- 
ness and power of the uc technique for this purpose. In 
the next section we briefly review this issue in which we 
combine the C-LSQ algorithm^ with the MC technique. 
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III. RESOLVING DISTRIBUTION OF 
RELAXATION TIMES 

A. Numerical procedure 

The complex dielectric permittivity, e = e' — le" , of 
materials in the distribution of relaxation times approach 

is expressed as follows 

e(a;) = e'{LO > lo^) (1) 

/oo 
g (log r) [1 + twr] ~ ^ d log T 
-oo 

where u is the angular frequency (w — 2ttv) of data in the 
frequency window with a lower bound ui and with a high 
bound u)-[ {oji <co<co-\)- r is the relaxation time of the 
process with the total dielectric strenght Ae. e'{LO > u>i) 
is the instantaneous permittivity value that contains all 
relaxations faster than the considered frequency window. 
Our numerical procedure extracts the unknown distribu- 
tion function g(r)^^ in Eq. (1) from the dielectric spec- 
trum without any a-priori physical assumption. The re- 
laxation distribution function must satisiy the condition 

/oo 
g(l0gT)dl0gT = 1 (2) 
-oo 

Before continuing any further, attention should be paid 
to the frequency dependent permittivity e{Lu) in Eq. (1) 
which involves an integral dependent on dlogr. In nu- 
meric analysis, this is undefined (or ill-mannered) be- 
cause one has a finite number of data points A'^ - in this 
case we have to use discrete r values such that the num- 
ber of r's (A^o) should not exceed N. So Eq. (1) is instead 
written as a sum of individual relaxations. 

£(W) = s'iiV > Wt) + Ef° A£i[l + lUTD,]-' (3) 

In this approximation, the Maxwell superposition prin- 
ciple with non-interacting independent dipoles (assump- 
tion of Debye ''^) with relaxation times td and dielectric 
strengths Ae are assumed. Then using CNLSQ, one can 
find pairs of [td, Ae]. However, the pairs would not be 
unique-we get different [td, Ae]-sets for different initial 
guesses or in other words the solutions are not orthogo- 
nal to each other. Our innovation in the method comes 
here, we introduce a pre-distribution of relaxation times 
(t). Fixing the Nq number of T/j-values, the optimization 
problem becomes linear and can be solved by C-LSQ-the 
only unknowns are Aej-values in this case which satisfy 

Ae = ^Ae, (4) 

Simultaneously, the values of Ae^ should not be nega- 
tive due to the physical insignificance of the results; we 

neglect any possible resonances in the dielectric spec- 
trum. This leads to assigning constraints on Aei, and 



the outcome of the least-squares for the Ae^ are all 
positive. Since there is a one-to-one relation between 
the pre-distributed td^ -values and the resulting dielec- 
tric strengths Ae^, the probability densities of relaxation 
times T£). are 

Sim.) = ^ (5) 

This is the resulting relaxation time spectrum, g(TD.), 
and Ae is defined in Eq. (4). However, keep in mind 
that a single optimization run for gir) does not mean 
anything-the solution is not unique and the time-axis is 
discrete. To overcome this, a MC procedure is introduced 
to obtain the [td , Ae] pairs-the optimization repeated 
many times (n)-for randomly selected td values^* us- 
ing a pre-distribution function. For the calculations the 
time-space can be assumed to be continuous. The distri- 
bution of the generated Tjji and the obtained Ae^ from 
C-LSQ optimization determine the relaxation time spec- 
trum, g(logrD). 

The obtained distributions are analyzed by means of 
comparing them with known distributions. We apply the 
Levy statistics^^'^^"^^ , which is used for interacting sys- 
tems in different research fields^*'^^"®^. The Levy stable 
distribution is a natural generalization (approximation) 
of the normal (Gaussian), Cauchy or Lorenz and Gamma 
distributions. It is used when analyzing sums of indepen- 
dent identically distributed random variables by a diverg- 
ing variance. Its characteristic function is expressed as^^ 

L(a;;AM,7,C)=^|exp{-|C(a;-M)r}l (6) 

Here, < 7 < 2. We need only four parame- 
ters; 7 characteristic exponent, /i localization param- 
eter, ( scale parameter and A amplitude. The spe- 
cial forms of Eq. (6) are the Gaussian [L(a:;; A, /i, 2, ()], 
the Lorentz or Cauchy [L(.t; A, /i, 1, ()] and Gamma 
[\-{x; A, ji, 1/2, Q] distributions. This equation could be 
used to characterize the properties of the obtained drt 
with A corresponding to the amplitute of the mean re- 
laxation time-its dielectric strength. 

B. Illustrative examples 

1. Delta function (distribution) 

A distribution in the form of a delta function yields the 
integral in Eq. (1) to be equal to the Debye relaxation 
function, g(logT) = 5i\ogT — logro) with the relaxation 
time tq. We have added 1% Gaussian random error to 
the data. The data is generated by assuming a relaxation 
with T = (27r)^^ s with Ae = 1. In Figure 3, the results 
from the dielectric analyses are shown. In the analysis 32 
relaxations are used and the Monte Carlo procedure is re- 
peated for 6400 times. The solid ( ) lines in Figure 3a 

are the fits obtained after the reconstruction of the orig- 
inal data. The drt is presented in Figure 3b is analyzed 
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using Eq. (6) . The characteristic exponent yields the dis- 
tribution is one of the special cases of Eq. (6), Lorentz 
or Cauchy distribution, 7 ~ 1 and ^ ~ 0. The Cauchy 
and the Levy distributions are illustrated with the solid 

( ) and dashed hnes ( ), respectively. Since the 

number of unknown parameters in the Cauchy is lower 
than the Levy's, the Cauchy distribution is used in the 
following discussions regarding the PPG data. Moreover, 
the Cauchy distribution is also one of the delta sequences 
C(r;AT,n)6780^ 

C{T;A,f,n) =i L(x;A ^, 1, C) (7) 
= ^n7r"^{l-|- [n(logr - logf)]^}"^ 

where A and n are fitting parameters and f is the most 
probable relaxation time. The actual relaxation time is 
also displayed at — logr = log27r = 0.8 as a vertical line 
in the figure. Eq. (7) can also be converted to the Breit- 
Wigner (also known as the Lorentz) distribution, which 
is a general form originally introduced to describe the 
cross-section of resonant nuclear scattering, by substitut- 
ing n with 2/ At, where Ar is the width at half maximum 
(FWHM). Although f and At are related to each other 
with Heisenberg's uncentainity principle (f At = /i/(27r); 
h is Plack constant), we do not consider any quantum 
mechanical effects during the generation of the data, and 
there is no such relation between two parameters f and 
At. For various simulations of the Debye relaxation 
yields 6.5 < n < 10 or 0.2 < At < 0.3, the range is 
strongly dependent on the number of the Monte Carlo 
loops, and the actual position of the relaxation in the 
experimental window. We have also used Eq. (6) to fit 
the results, this equation was more successful for lower 
values of the drt, at high values there is no difference 
between the general Levy and the Cauchy distribution. 
Summing up the obtained distribution is verified to be a 
delta sequence, C(t; A, f, n) = S{t — td) 

The comparison between the generated data and the 
extracted distribution illustrates a line broadening^^, 
which is due to several factors. First of all we have added 
some Gaussian noise to the data, although this is not as 
significant as the other factors, it infiuences the distri- 
bution. Secondly, the td values are selected randomly 
from a log- linear distribution [see Ref.^ for details]. If a 
biased distribution is used i. e. based on the derivative of 
e' with respect to frequency (dloge/dlogw), the width 
of the distribution get narrower. However, it is nontrivial 
to find the initial distribution. If dloge/dlogw is used 
the resulting distribution is not smooth enough due to 
numerical derivation problems. The numerical integra- 
tion in Eq. (3) is the other reason, which because of the 
calculation-speed considerations is based on a simple nu- 
merical summation with unit base of each m generated. 
It can be improved by altering the integration routine. 
Finally, we can further improve the solution by increas- 
ing the MC loops, which is a statistical improvement. As 
a final remark, one can also argue the significance of a 
delta function in the concept of theory of distributions 
and real systems in the Nature^^. 
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FIG. 3: The Monte Carlo analysis of Dobyc relaxation with 
relaxation time t = (27r)~^ s (— log(T) ~ 0.798) and Ae = 1. 
(a) The real and imaginary parts of the dielectric susceptibil- 
ity, X = £ — e' {ijj > ijj^') and (b) the DRTcalculated with the 

presented procedure. The smooth solid ( ) line in (b) is 

the fit obtained from the Lorentz (Cauchy) distribution with 
n = 8.165 (FWHM = 0.245), 'logf = -0.801 and A = 0.018. 

The dashed line ( ) in (b) is the Levy stable distribution, 

L(logT; 3.768, -0.801, 0.834, 8.167). 

2. Box distribution 

Now, we apply the method presented to a box dis- 
tribution. The box empirical distributions has been de- 
scribed in detail for mechanical relaxations^^'^* and in 

the Frochlich molecular theory of relaxation^''. The box 
distribution is defined by the following equations 

H{logT) = constant for ti < t < T2 

if(logT)=0 for r>T2;T<Ti (8) 

The results are presented in Figure 4a and 4b, which 
show the reconstructed and obtained drt, respectively. 
The original box distribution is also presented in Fig- 
ure 4b as thick rectangle which indicate a good aggree- 
ment with the DRT. We have also applied Eq. (6) to 
fit the DRT. The result is shown as a dashed line ( 
), L(logT; 0.553, -0.788, 4.740, 1.019). The generalized 
spectral line expression is successful to model even broad 
spectral lines as the box distribution. 

IV. EXPERIMENTAL 

PPG of weight-average molecular weight 4000 (pur- 
chased from Polysciences Inc., USA) was used after high 
vacuum drying at 10 |xPa at 80 °C for 24 h. The cell was 
mounted in an inert atmosphere. The dielectric measure- 
ments on PPG samples were performed at temperatures 
between 248 K and 344 K in an angular frequency co range 
27rl0^ s-i - 3.67rl0^ s"\ using a Hewlett Packard HP 
4291 A HF impedance/material analyser together with 
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FIG. 4: The Monte Carlo analysis of the box distribution, 
(a) The real and imaginary parts of the dielectric suscepti- 
bility, X = e — e' (uj > uj-\) and (b) the DRT calculated with 

the presented procedure. The thick solid ( ) line in (b) 

is the original box distribution. The dashed line ( ) is 

L(logr; 0.553, -0.788, 4.740, 1.019). 



a Novocontrol BDS 2100 sample cell and a Novocontrol 
2200 RF extension line. The amplitude of the signal was 
0.1 V. 
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V. APPLICATION OF THE METHOD 

Parameters used in the MC analyzing procedure were 
as follows 

iV > 64 ATo = 32 n = 10000 

The number of td -values were 3.2 x 10^ in the calcu- 
lations. Computations were performed on a Pentium 4 
(2.4 GHz LinuxPc) with 1 Gbyte memory. Analyzing 
each dielectric data set takes between 3 to 4 h with the 
parameters selected above. 

It is very important to assign an instantaneous per- 
mittivity value to start the analyses-relaxations over w\ 
arc summed in this permittivity value, here oj-^ is the 
highest angular frequency used in the experiments. This 
is not actually just a constant number substracted from 
the real part of the dielectric permittivity since the data 
(the dielectric susceptibility) must afterwards satisfy the 
Kramcr-Kronig relation. For this reason the procedure in 
Str0mme-Mattsson et al. ^ has first been adapted, how- 
ever, in some data sets this procedure has not produced 
good-fits at high frequencies. Therefore, a new preprocess- 
ing method was introduced in which the high frequency 
part of data was curve-fitted by a Havriliak-Negami em- 
pirical formula [Eq. (9)]"^^ with a CNLSQ algorithm with 
a fast relaxation time r « Later, the data-set is 
extrapolated by this empirical formula, which signifi- 
cantly improves the DRT-spcctra. Obtained drt is an- 
alyzed by Eq. (7) and the results are compared with 



FIG. 5: (a)Roal e' and (b) imaginary e" parts of permittivity 
versus logarithm frequency v at various temperatures. 

those obtained from the curve-fitting analyses using the 
Havriliak-Negami expression, 

e = e'{uj>iJ^) + Ae[l + {iujT)°']-l^ (9) 

where, Ae, r, a and /? are the fitting parameters. 

VI. RESULTS AND DISCUSSIONS 

Figure 5 shows the real and imaginary parts of the di- 
electric permittivity for PPG over the available frequency 
range. Each curve corresponds to a different tempera- 
ture. The a-relaxation is shown as a distinct peak, and 
the onset the /3-relaxation can be observed at the low- 
est temperatures. At temperatures over 300 K the two 
peaks overlap, producing an increase in the magnitude of 
the dielectric strength. At temperatures above 300 K a 
weak a'-relaxation may be seen at low frequencies. The 
frequency of the dominant a-peak increases with temper- 
atTire and its amplitude slightly decreases with increas- 
ing temperature. This behavior is characteristics of many 
molecular systems, i.e., sec for example. Moon ct al.^^ 
for PPG and Dixon et al. Schonhals et al. ''^ and Angell 
et al. for other glass-forming liquids. Although curve- 
fitting approaches using a Fourier transform of streched 
exponentiaF^'^^ or Havrialiak-Negami"^^ as well as its 
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FIG. 6: Dielectric permittivity at high frequencies e'{u) > lo^) 
(•), double maximum dielectric loss 2eJ^iax (■) and As (x) as 
a function of temperature. Lines are drawn to guide the eye. 



derivatives would be an appropriate method for the anal- 
ysis, it does not produce valuable information regarding 
individual molecular processes other than the most prob- 
able one^^. In addition, it is worth mentioning that at 
temperatures higher than the glass-transition tempera- 
ture, it is not easy to observe the /3-relaxation in some 
molecular systems since it is hindered by the stronger a- 
rclaxation in the crossover region. Finally, the dielectric 
loss of PPG does not exhibit any sign of ohmic-losses at 
high temperatures and at low frequencies. 

Figure 6 shows the instantaneous dielectric permittiv- 
ity as a function of absolute temperature at high fre- 
quencies e'{iv > w-f), double maximum loss and 
calculated dielectric strength Ae. The permittivity value 
fluctuates between 1 and 1.05 at temperatures lower than 
300 K. It indicates that the sizes as well as the number 
of dipoles with fast relaxation times are small at frequen- 
cies higher than 1 GHz. At higher temperatures instan- 
taneous permittivity increases-this means that a couple 
of relaxations are finalized and they contribute to in- 
stantaneous polarization. The estimated total dielectric 
strength Ae, on the other hand, neglecting the first data 
point (the first data point could be due to virgin sample, 
which experiences the temperature and electric field the 
first time), decreases as the temperature increases. Tak- 
ing into account the Debye relaxation model^^'®^, which 
states that Ae = the discrepancy between the 

calculated Ae and double the measured maximum loss- 
peak 2 • e^jjx proves that the relaxations at temperature 
lower than 320 K are broad (non-exponential decay func- 
tion) and can be modeled by the drt. The data also in- 
dicates that there can be several single relaxations with 
relaxation times close to each other. However, at tem- 
peratures higher than 320 K these two values get close 
to each other indicating that the relaxation is nearly- 
exponential. The value of 2 • e[[-^^^ must not exceed the 
dielectric strength, but at the highest temperatures this 
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FIG. 7: (a) Dielectric susceptibilities x versus angular fre- 
quency u} and (b) the product of dielectric strength and drt 
Ae X g(logr) obtained for temperatures 322, 333, and 344 K. 
The symbols (o) and (o) are the real and imaginary parts of 
the susceptibilities after subtraction of £'{u! > W|)-they are 
the measured data from experiments. Solid ( ) lines rep- 
resent the reconstructed resposes using the distribution func- 
tions g(log r) for the real and imaginary part of the dielectric 
susceptibilities, respectively. 



condition is not satisfied, which is due to lack of enough 
data points at high £requencies-we have no [r, Ae] pairs 
at high frequencies from the calculations. 

The analyses by means of the drt for three tem- 
peratures, 322, 333, and 344 K, are displayed in Fig- 
ure 7. In the figure, curves in the upper graphs (Fig- 
ure 7a) are the real and imaginary parts of the suscep- 
tibilities, x(ci;) = e(a;) — e'{uj > a)|), and the curves in 
the lower (Figure 7b) are the logarithm of the product 
of the total dielectric strength Ae and the obtained drt 
g(logT). The reconstructed dielectric susceptibilities are 

illustrated as solid ( ) lines for real and imaginary 

parts. The symbols in the figure are the measured sus- 
ceptibilities X after substruction of instantaneous permit- 
tivity e'{co > u>-[). The drt is presented as histograms. 
Focusing on the three peaks obtained from the Monte 
Carlo analysis in the drt spectra. Figure 7b, the most 
probable relaxation time behavior can be seen as stated 
by Angcll et al.^''. Neglecting the weak relaxation at 
logr ~ —7.5, and just considering the data in the an- 
gular frequency region 100 Ms~^ > w > 20 Gs~^ - one 
could assign a relaxation time using a curve-fitting pro- 
cedure which would probably yield only one relaxation 
time (the superposition of the relaxations as observed 
in Figure 7b). As a remark, our analyses indicate that 
there exist a couple of relaxations taking place in PPG. 
Interestingly only one of the peaks move with tempera- 
ture. The peaks at 10 Gs"""^ are due to the preprocessing 
of the data since the spectra should start with a Debye 
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FIG. 8: The product of dielectric strength and the drt 
[Ae X g(r)] versus logarithm of relaxation time r at differ- 
ent temperatures. 

relaxation^. However, this fast relaxation, which shows 
no temperature dependence, could also be treated as a 
boson peak, which seems to be a general phenomenon of 
disordered molecular systems and not connected with a 
glass transition^ ^. 

The time-temperature superposition states that by 
heating the material, the dielectric loss peak moves to 
higher frequencies. This movement is related to the 
thermally activated relaxation processes, that makes it 
possible to estimate the activation energy for the par- 
ticular loss. In Figure 8 DRT versus relaxation time 
logr is plotted at different temperatures. The relax- 
ations at very short times (logr ^ —10) are due to the 
experimental conditions since there exists a first relax- 
ation process in the material with the application of the 
voltage^. This relaxation could also be confused with the 
boson peak, which shows no temperature dependence as 
stated by Donth-"-^ and Dionisio et al.^^. The amplitutc 
of the high frequency relaxation increases as the tem- 
perature is increased due to the movement of the di- 
electric data to higher frequencies where they show a 
frequency-temperature superposition and interfere with 
a third peak not sensible to temperature. Besides, two 
peaks at low temperatures becomes one peak after 300 K. 
The two relaxations at low temperatures with slow re- 
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FIG. 9: Relaxation-image for poly (propylene glycol). The 

solid lines ( ) represent the fitted curves of Ref. 71 for a- 

and /3-relaxations, and of Ref. 13 for the a'-relaxation, respec- 
tively. The symbols (•) represent the segmental relaxation for 
poly (propylene oxide) taken from Ref. 8. 

laxation times (logr ~ [—8.5 ... — 7]), do not mix with 
each other on the crossover region^^, which is clearly ob- 
served in the figure. The method applied to resolve the 
DRT gives a new insight into this region in glass-forming 
liquids and polymers. The amplitude of the slowest re- 
laxation decreases just before the crossover temperature 
Tc (approx. 275 K), however, the amplitude of the inter- 
mediate relaxation increases its amplitude as the temper- 
ature increases. The position of this peak starts to move 
to faster times as temperature values are over 275 K. 

In order to illustrate the temperature dependence of 
dielectric relaxation for PPG, a different representation 
is shown in Figure 9. In the image, the dark spots are 
the hills of the relaxations and white regions are plateaus 
where no relaxations are observed. The lines in the fig- 
ure are taken from Arrhenius plots of previously pub- 
lished papers^'^'^^. The solid line ( ) in the figure 

represents the fitted-curve to a VTF equation adapted 
from Moon ct al. for the a-relaxation. Our tempera- 
ture interval is just above the temperature region that 
Moon et al. "^^ had in their experiments. At low tempera- 
tures the obtained drt coincide with the fitted curve in- 
dicating that the main diagonal peak is the Q-relaxation. 

The solid line ( ) in the figure is an Arrhenius fit to 

the /3-relaxation again from Moon et al. As stated 
by Moon et al.^^, near 240 K the a- and /3-relaxations 
separate, however, the drt obtained demostrate that 
this separation or merging of the two relaxations occur 
around 275 K. The discrepancy could be due to the 
sample preparation and the position of the glass transi- 
tion temperature which affects the cooperative processes 
in the sample. Moreover, the separation region in the 
relaxation-image is broad and the drt analysis probes 
this region better than the conventional methods. The 
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FIG. 10: DRT at 299 K. The data obtained from the presented 

analysis is displayed with error bars. The dotted lines ( ) 

arc the three Lorcntz distributions fitted to the fast, interme- 
diate and slow relaxations. The solid line ( ) shows the 

sum of these Lorentz distributions. 



proximity of the /3-relaxation peak induces a change in 
the a-relaxation, which is clearly observed as a repulsion 
of two peaks in Figure 8. It is interesting to sec this ef- 
fect as if the /3-relaxation slows down the a-relaxation 
as they approach each other around 275 K. In addition 
the amplitutcs of the two relaxations behave opposite to 
each such that the amplitude of the /3-rclaxation increases 
while the amplitude of the a-relaxation decreases as they 
approach to each other. This illustrates that the coop- 
erative motion of the segments loose their character and 
some part of them join the local movement of the struc- 
tural units, however, the activation energy of this new 
process is higher than the one for the local movements 
(/3-relaxation) . 

We should always keep in mind that the VTF equation 
uses the most probable relaxation time. At high temper- 
atures, the VTF equation is therefore not valid which is 
also expected after the merging of two relaxations. Al- 
though there is a good agreement between the data of 
previous investigations and ours, the peaks (dark spots) 
of the DRT is not coinciding exactly. The diagonal move- 
ment of the dominant peak in Figure 8 is significant and 
there are two horizontal spots that contribute to the diag- 
onal main sequence. The data of Mijovic et al. ^ are also 
illustrated in the figure with symbols (•), which represent 
the segmental mode relaxation for for poly (propylene ox- 
ide), and it coincides with the a relaxation of ppg as 
reported by Moon et al. . 

The obtained DRT as presented in Figure 8, are ana- 
lyzed by the Lorentz or the Cauchy distribution [Eq. (7)] 
to find the more exact positions and widths of the in- 
dividual DRT, Ae X g(r) = C(r;^,-f,n). An example 
of this curve-fitting procedure is displayed in Figure 10, 
where the DRT results are shown with errorbars and the 
fitted Lorentz distributions are the solid lines and Table 



TABLE I: The Cauchy or Lorentz distribution [Eq. (7)] pa- 
rameters for the DRT at different temperatures. 
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I presents the parameters of the distributions at differ- 
ent temperatures. In the table relaxations in each tem- 
perature arc separated into three different contributions, 
fast, intermediate and slow, respectively (also presented 
in Figure 10). Similarly the data are fitted by Eq. (9) and 
the fitting parameters are given in Table II. As men- 
tioned previously, while discussing Figure 6, the curve 
fitting procedure with the Havriliak-Negami expression 
confirms that the dominant relaxation becomes Debye- 
like as the temperature approaches to 340 K. However, 
this procedure is neither probing into the crossover region 
or resolving the present two relaxations. In addition the 
Havriliak-Negami curve fitting procedure yields permit- 
tivity values at high frequencies e'(a; > Wj) lower than 
one which is not physical for PPG. Returning back to 
Table I, it is illustrative to see the change in the ampli- 
tude of the fast relaxation as the temperature is increased 
in the experiments. While comparing the amplitudes of 
the intermediate and slow relaxations for temperatures 
around the crossover temperature, we can decompose 
the overlapping dynamics of a and /? processes, using 
the probability theory. The presented analysis technique 
can provide a deeper insight of processes taking place in 
molecular systems such as dipole-dipole interaction. 

Last but not least we display in Figure 11 the Ar- 
rhenius plot of PPG. In the figure the logarithm of in- 
verse relaxation times (t~^ = w) are presented. The 
results obtained from both the Havriliak-Negami empir- 
ical formula and the DRT method are used. The open 
symbols (o, □ and o) show the DRT analyses (Table I). 
The filled symbols (★) are the relaxation times obtained 
from the curve- fitting of the empirical formula (Table II). 
To give some more insight, data from the literature are 

also added; the dashed ( ) and chain (- • -) lines 

are from Moon et al. and the dotted line ( ) 

is from Schonhals and Kremer^^. The graph is divided 
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TABLE II: Fitting parameters obtained by applying the 
Havriliak-Negami empirical formula in Eq. (9). 
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into four sectors, which are used in the discussion below. 
The origin of the sectors is the crossover temperature 
and the corresponding relaxation rate (logWc)- In sec- 
tor 1, there is only the trace of the j3 relaxation and 
the DRT analysis yields relaxation rates little higher than 
stated in the literature^^ . In sector 2, there is a couple of 
points that should be mentioned. First, the curve-fitting 
method and DRT result in different relaxation time con- 
stants. The DRT result in a linear Arrhenius relation, 
logw = a + h/{T) with a = 13.47 and h = -1374, while 
the curve-fitting method show a good agreement with 
the VTF equation reported in literature''^. One should 
keep in mind that our experimental interval is just out- 
side of those in Moon et al.''^, and the line drawn is an 
extrapolation of their results. In sector 3 there is only 
the weak slow relaxation, which is not as significant as 
the others in the DRT spectra. The relaxation could be 
associated with the a'-relaxation. However, we should 
keep in mind that the relaxation has a weak amplitude 
and data is just at the border of our experimental win- 
dow. Lastly in sector 4 the a relaxation is present and 
it vanishes as it approaches the crossover temperature 
Tc- Although the Havriliak-Negami empirical formula 
produce valuable results, it is not sufficient to give in- 
sight into the crossover region as the DRT analysis. It 
should be noted that the behavior of three relaxations 
shown in Figure 11 are very similar to those illustrated 
for polyethyl methacrylate in page 29 of Donth^'' also in 
Ref. 11, for poly(alkyl methacrylate)s of Beiner-"^^ and 
for other glass-forming systems^^'^^'^^'''''. 



VII. CONCLUSIONS 

A numerical method is applied to dielectric data of 
PPG to extruct the DRT. To demonstrate the potential 
of the method, the drt of the known distributions, delta 
function and box, were obtained from generated data. 
Later the same procedure was applied to experimental 
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FIG. 11: Arrhenius plot of PPG. Symbols show the result of 
the present study. Open symbols [(o)-fast relaxation, (□)- 
intermediate relaxation, (o) slow relaxation] are relaxation 
rates obtained from drt analyses (Table I) and filled symbols 
(★) arc the curve fitting results of the employed Havriliak- 
Negami empirical formula (Table II). The points represented 
by (-I-) symbols are the relaxations with large amplitudes in 

Table I (^ > 5 x 10"^). The sohd ( ) fine is a straight 

line fitted to the intermediate relaxation at high temperatures 

(T > 275 K). The dashed ( ) and chain (- ■ -) lines are 

the VTF and Arrhenius fits from literature^^ for the a and /3 

relaxations, respectively. The dotted ( ) line is the VTF 

fit for the a'-relaxation^^ . The figure is divided into four 
sectors taking the crossover point [10^ /Tc, logWc] as center. 

data of PPG. The method was able to find a unique DRT 
for given dielectric data. It was illustrated that for molec- 
ular systems we were not limited to analyze average or 
most probable relax;ation times. However, new questions 
arise for the origin of the obtained drt. 

PPG showed two significant DRT which are visible be- 
tween 260 — 277 K. At the higher temperatures there was 
only one visible peak in the DRT spectra. The most sig- 
nificant relaxation below 275 K was associated with the a 
relaxation. The second significant relaxation on the same 
temperature region was associated to the /3- relaxation. It 
was observed that the a- and /3-relaxation separation (as 
well as merging) (crossover region) was around 275 K. As 
the temperature increased above the crossover tempera- 
ture, only one relaxation was significant, and it had a lin- 
ear Arrhenius relation. It appeared as if the /^-relaxation 
slowed-down the a-relaxation which proved the strength 
of the applied analysis method over the traditional curve- 
fitting procedures. Comparison of the DRT approach to 
those of the curve-fitting one have especially pointed out 
the potential of the DRT technique in the case of the 
crossover region. The process at temperatures above the 
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crossover temperature was not actually a continuation 
of the a relaxation, which shows a VTF type tempera- 
ture dependence, but it was clearly a continuation of the 
/3-relaxation with a higher activation energy. The com- 
parison of the Arrhenius- or relaxation-image with data 
from literature verified that the method was successful 
for low temperatures, a- and /3-relaxations. Although, 
the fast relaxation resolved from the dielectric data coin- 
cides with the boson peak speculation in disorder molec- 
ular systems^^, other indications are needed to support 
the results. 

The DRT illustrated the importance of interaction be- 
tween different processes in molecular systems such that 
as the merging of the a- and /3-relaxations. Their inter- 
action could be observed by the presented method. The 
cooperative motion of the chain units (a process) creates 
space for the local motion (/3-process). As the sample 
is heated only the local motion of polymer units sur- 
vives as a detectable relaxation process, and it shows an 
Arrhenius temperature dependence. It is important to 
continue dielectric measurements on PPG samples with 
smaller temperature intervals and a broader frequency 
window to better understand the nature of molecular pro- 
cesses taking place. 

Finally, we in this article gave our attention to one 
of the ignored subjects, the drt in condensed-matter 
physics. It is illustrated through out the text that it can 
be used to obtain valuable information from the dielectric 
data, which is not possible otherwise by classical dielec- 
tric data analyses approaches whose main assumption is 



using the average or most probable relaxation time for 
a given dielectric loss peak. As a metaphor, it is impor- 
tant to state that in early days of astronomy well-known 
pioneers in the field used brightness of stars which were 
sum of all light from stellar-objects. Later researchers 
used the optical spectra to obtain abundances of differ- 
ent elenicints. which in return supplied valuable informa- 
tion for gaining a better understanding of stellar-objects. 
Well, traditional methods to analyze dielectric data only 
considers the shape of the dielectric data, and do not 
extract what is underneath. We, therefore, think that 
the applied numerical method acts in a similar way to 
the optical spectroscopy for astronomy, and resolves dif- 
ferent relaxations, "spectral-lines", hidden in the dielec- 
tric data. The next step is then to relate these individ- 
ual relaxations to structure and abundances of molecular 
species in the material, and better understand whether 
the lines are homogeneously or inhomogeneously broad- 
ened. 
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